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Collapse models predict the spontaneous collapse of the wave function, in order to avoid the 
emergence of macroscopic superpositions. In their mass-dependent formulation, they claim that the 
collapse of any system's wave function depends on its mass. Neutral K, D, B mesons are oscillating 
systems that are given by Nature as superposition of different mass eigenstates. Thus they are 
unique and interesting systems to look at, for analyzing the experimental implications of such 
models, so far in agreement with all known experiments. In this paper we derive — for the single 
mesons and bipartite entangled mesons — the effect of the mass-proportional CSL collapse model 
on the dynamics on neutral mesons, including the relativistic effects. We compare the theoretical 
prediction with experimental data from different accelerator facilities. 

PACS numbers: 03.65.-w, 03.65.Ud, 03.65.Ta 
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Flavored neutral mesons, those with a net non-zero strangeness, charm, or beauty content, are among the most 
fascinating systems in elementary particle physics. Particle and antiparticle are distinguished only by the flavor 
quantum number, exhibiting the phenomenon of flavor oscillations in their time evolution [l|, 0] (see also 
the following we will focus on K-mesons, but our conclusions can be easily generalized to the other flavored meson 
systems, as discussed at the end of the paper. 

Flavor oscillation is a direct consequence of the superposition principle of quantum mechanics. Nowadays many 
scientists question its validity [lT| - [l4| . and several experiments have been performed [l5l-f20j or proposed, which 




challenge it. On the theoretical side, models of spontaneous wave function collapse [2l|-|27[ explicitly predict that the 
superposition principle is valid only at the microscopic scale, while it gradually breaks down when moving towards the 
macroscopic scale. Therefore, by altering the standard quantum dynamics, collapse model predict a different behavior 
for kaon oscillations. Aim of this paper is to compute such an effect, and compare it with available experimental data. 
To this end, we will use one of the most popular collapse model, the mass-proportional CSL (Continuous Spontaneous 
Localization) model [H, [2^, and will perform the calculation of the probabilities, that are observed in experiments, 
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I. INTRODUCTION 
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to second order perturbation theory and from that deduce the higher orders. Comparison with experimental data will 
allow to conclude whether current or future experiments are sensitive to this effect. 

One of the difficulties in computing the effect of collapse models on kaon oscillation is that these models act on the 
spatial part of the wave function — since the collapse is supposed to localize wave functions in space — while oscillations 
occur for the internal degrees of freedom. What happens then, is that the noise responsible for the collapse of the 
wave function acts like a random medium through which the particles propagate, modifying their evolution and also 
the oscillatory behavior. By resorting to standard quantum field theoretical tools, we will compute analytically such 
an effect. 



II. KAON PHENOMENOLOGY AND THE MASS-PROPORTIONAL CSL MODEL 

We start by recalling why kaons oscillate. For convenience of the readers we keep the discussion here short and 
more details can be find in the Appendix A and references therein. The basic observation is that the mass eigenstate 
are different from the flavour eigenstates. For kaons, there are two mass eigenstates, the short state \Ks) and the 
long state \Kl), as well as two flavor eigenstates, \K^) and \K^)- In this paper, we work in the approximation that 
there is no CP violation: this implies that \Ks) and \Kl) are orthogonal^. The relation between the mass and the 
flavour eigenstates is given by: 

During the time evolution, the mass eigenstates change by acquiring different phase factors, depending on their mass^. 
The phenomenon of flavor oscillation arises because what we measure are not the mass eigenstates, but the flavor 
eigenstates, which are superpositions of mass eigenstates. The different phase factors in front of the mass eigenstates 
change the superposition, making it possible, for example, to start with a \K^) and end up with a j^'^). 

According to collapse models, the time evolution of the mass eigenstates is different compared to the one given by 
standard quantum mechanics. Hence we expect to see a different behavior in kaon oscillations. According to the mass 
proportional CSL model [l^ [2^ , the evolution of the state vector is given by the non- linear equation: 



-Hdt+ y-^ I dx (A/(x) - {M(x)))dWtU) ^ / dx (M(x) - {M(x))f dt 

h mo J 2mQ 



10*) , (2) 



where (M(x)) := {(pt \ (pt)- Here H is the standard quantum Hamiltonian of the system and the other two terms 

induce the collapse of the wave function in space. The mass mg is a reference mass, which is taken equal to that of a 
nucleon. The parameter 7 is a positive coupling constant which sets the strength of the collapse process, while Af (x) 
is a smeared mass density operator: 

M (x) = Y.m,N, (x) , N, (x) = j dyg (y - x) ^] (y) (y) , (3) 

where ^p'^ (y), ipj (y) are, respectively, the creator and annihilation operators of a particle of type j, namely having 
mass mj and spin s, in the space point y. Neutral kaons arc spin zero particles, thus the spin will be of no relevance 
for the following calculations. The smearing function g{x.) is usually taken to be a Gaussian: 

ff(x) = , J , e--'/^^i , (4) 



27rrc)" 

where rc is the second new phenomenological constant of the model. The standard numerical value of this correlation 
length rc is [Hill ill: 

rc = 10~^cm, (5) 



^ This is justified by the smallness of the CP violation effect in neutral kaons, which gives rise to a small (of the order of 10"'^) odd/even 

CP impurity in the Kg/Kj^ states, and to a small non-orthogonality between them. 
^ Indeed, this is not exactly true: if we take into account that kaons decay in time, there is also an exponential dumping factor, besides 

the phase factor. We will consider this property later. Here we want to explain the basic idea behind the oscillation phenomenon. 
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while, in the hterature, two different values for the collapse strength 7 have been proposed. The first value has been 
originally proposed by Ghirardi, Pearle and Rimini 

7 = 10-3°cm^s-i (6) 

in analogy with the GRW model [2l|. The second value has been proposed by Adler, inspired by the analysis of the 
process of latent image formation according to collapse models, and amounts to (25j 

7 = lO^^^cm^s-^ (7) 

Finally, Wt (x) is an ensemble of independent Wiener processes, one for each point in space. 

Working with non-linear equations is notoriously difficult. As shown e.g. in Ref. [231 the experimentally testable 
predictions of the model — when averaged over the noise — do not change if the real noise Wt (x) is replaced by an 
imaginary noise iWt (x). In this way, one loses the collapse properties, i.e. the non-linear terms, of the equation. 
However, the advantage of having an imaginary noise is that the evolution is described by a standard Schrodinger 
equation with a random Hamiltonian: 

iJ.oT = H-h^Y. — I dyw{y,t)^]{y)i;j{y) (8) 
^ mo J 

3 

where 

w{y,t) = j dxg(y-x)et(x), (9) 

and Ct(x) = dWt(x)/dt is a white noise field, with correlation function E[^t(x)^s(y)] = 5{t — s)5{x. — y). As such, 
w(x,t) is a Gaussian noise field, with zero mean and the correlation function: 

E[w{^,t)w{y,s)] = ,5(t-,s)F(x-y), F(x) = -^^—e-^'/^^l. (10) 



In the following, we will be a bit more general without making the calculation more complicated, and we will assume 
that the noise has a general time correlation function /(s) instead of a white noise time-correlator S{s). 

In our analysis, we treat the kaons as non-relativistic particles, which is in accordance with the experimental 
situation at some acceleration facilities. Accordingly, the Hamiltonian for the mass eigenstates is given by (j = S,L] 
short or long) 

H{t) = j d^n.ix) (11) 



j=S,L 



with 



^2 

Hj {x) = m^cV] (x) V'j (a;) + ^^Vt/.] (x) • {x) - h^J^w {x) ^] {x) V'j (x) (12) 



J3 



and 



: Af^ {X) 
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7^,-7^) . (13) 

The Hamiltonian includes two contributions: the standard Schrodinger term T-Ls {x) = T-Lg (x) -f- T-Lg {x) and the term 
M {x) = {x) +Af^ (x) which accounts for the collapse. Note that in Hs (x) we have included also the mass energy 
terms. These terms are usually ignored, since they lead only to a constant shift of the energy, which implies no 
observable effect. In the specific situation of kaons, where we have a superposition of two different mass eigenstates, 
it is fundamental to keep it. 

Indeed, for most investigations of the kaon phenomenology the kinetic part is not relevant and one introduces an 
effective Hamiltonian with two different mass eigenstates, i.e. the system is treated as a two state system analogously 
to spin-i particles [1, Q and the two relevant bases — if CV violation is neglected — are the strangeness (flavor) basis 
and mass eigenstate basis. For convenience of the readers we provide a summary in the Appendix A. 
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Since the effect of the CSL noise is to localize the wavefunction in space, we cannot simplify our model by neglecting 
the kinetic part and we need to work in a standard field theoretical background. In order to compute the effect of 
collapse models on their evolution, we must specify therefore the spatial part of the wave function, the form of the 
wave function in space, which will be localized via the background field. This means that, as said before, the change 
in the oscillatory behavior of kaons, according to the collapse dynamics, is an indirect, not a direct, effect of the 
collapse process. 

In the following, we focus our attention on computing the probability that, starting from a ji^*^) state at time t ~ 0, 
we end up in a \K^) state at a later time. The others possible transition probabilities (jii'') — > l^"); |^") — ^ 1^") 
and ji^*^) — > 1^*^)) '^^^ t)C computed in a very similar way. In order to keep the computation as simple as possible, 
we assume that the initial state is a plane wave with definite momentum p^. Therefore, the quantity we wish to 
compute is: 

Pko^Ro {P^) ^Y.^\{K",Pf |C/(i)|i^°,P.)|' (14) 
Pf 

where we assumed that also the final state is a momentum eigenstate, and we sum over all the possible final states. 
Here, E denotes the stochastic average with respect to the noise of the background field. It is convenient to express the 
matrix elements in terms of the mass eigenstates, instead of the strangeness eigenstates, since they are the diagonal 
states of the Hamiltonian and, therefore, provide the simplest form of the time evolution: 

{K°, Pf \U (t)| p,) = J2 {K., Pf \U it)\ K„p,) = J2 {Kj,Pf \Uj {t)\ K„p,) 

i-j 3 

EE Y.ajl5*T,{pf,p,,t) , (15) 

j 

where i,j ~ S,L and as = ctL = I3l = l/"\/2, Ps = — 1/\/2 [The parameters a, /3 change accordingly for the 
other transition probabilities] . The second equation is due to the fact that since the structure of the Hamiltonian is 
given by H (t) — Hs (t) + (t), namely not mixing of the mass eigenstates, the time evolution operator factorizes 
U (t) = Us (t) ^ Ul (t). This implies that if a mass eigenstate is produced at time t = it persists its identity for any 
later time point. This is strongly supported by experimental data. 
Substituting Eq. ([15]) into Eq. ([14]) we get 



Pro^Ro iP^, i) - E «^-/3;«fe/3feE ^ K (P/, P» t) T, (P/, P., t)]. (16) 

= Pkj ip^;t) 

We now compute this probabilities. 

III. DERIVATION OF THE PROBABILITIES FOR ONE-PARTICLE STATES 

We cannot solve the equations of motion exactly, because of the noise term in the Hamiltonian. Therefore, we 
move to the interaction picture and apply a standard perturbative approach. We treat {x) as the unperturbed 
Hamiltonian, and Af (x) as a perturbation. This is certainly a very reasonable assumption since the noise coupling 
constant 7 is very small. 

A. The Interaction Picture 

In the interaction picture, the states evolve as follows [l^ 

\i;t)j = Ui{t,0)\i;o)j , (17) 

with 

Ui {t, 0) = 1 + i^—j dh dt2... dUNi ih) Ni (h) ...Ni it„) (18) 
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given via the well known Dyson series. The relation between the evolution operator in the Schrodinger picture and 
the corresponding one in the interaction picture is given by 

U{t,0) = e-^"^*Ui{t,0). (19) 

The fields, on the other hand, evolve accordingly to the free Hamiltonian, therefore we expand the wave function into 
a superposition of plane waves 

^,(x) = ^^6ke-^(^'='-''--) . (20) 

In order to avoid possible divergences coming from the fact that we are working with plane waves, we quantize our 
system in a box with side L by imposing periodic boundary conditions. This implies that the components of the 
wave vector k take only discrete values kj — ^jj^n-kj, with Ukj integer. The energy is ~ mc^ + ^ and is the 
annihilation operator of a particle with momentum k. Since kaons are bosons, their operators have to satisfy the 
commutation relations [6k,&j(./] = <5k.k'- The total Hamiltonian can then be written as 



H = J2Ekblbk. (21) 



k 

At the end of the calculation, we take the limit L — > oo, which amounts to make the substitutions 

+00 



— — OO ^ 



r 

k— — OO 



We have now introduced all the necessary elements for computing the matrix elements and the transition probabilities. 



B. Computation of the Transition Amplitudes 

We start by focusing our attention on the computation of the transition amplitude of a certain mass eigenstate 
Tj {pf,pi,t). Since this computation is independent of the mass eigenstate j = S, L, we drop this index in this section. 
The interesting expression to compute is 



T ipf, p„ i) = {pf\U it) |p,) = {pf\e-^"^*Ui (t, 0) |p,) = e-i'^f'inibp, Ui (t, 0) 



(23) 



where Ef = ^ and \n) is the vacuum state. The perturbative scheme goes as follows. We keep terms up to second 
order, which — when averaged over the noise — give the first order contribution to the oscillation. The reason is that 
the corrections are proportional to the average of products of the noise computed in different space time points, which 
are not zero only when the noise appears at least twice in the matrix elements. Therefore we have 



Ui (t, 0) ~ 1 



Jo 



1 



ti 



dtiNi (ti 

Accordingly, the transition probability becomes: 

T (p/, p„ t) ~ e-i^f' [rW (p/; p.; t) + T^^^ (p/; p,; t) + T^^) (p^; p,; t) 



(24) 



(25) 



where each term corresponds to one of the first three terms of the Dyson series. It is possible to give a representation 
by means of Feynman diagrams for each one of these three terms: 



where the solid lines refer to the kaon and the dotted lines to the noise. The first term, T'^*^-' , is trivial: 

TW(p/;p,;t)^(l]|6p,6^Jf7)=5p,,p,. (26) 
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The second term is given by: 

T^^HPf■,P^■,t) = ^Vl^ I dh I d^iw{xi){n\bp^4{xi)^iixi)bljn). 



(27) 



Using the plane-wave expansion of the fields given by Eq. (j20p . and a similar expansion for the adjoint of the field, 
we get: 

{n\bp,4 (xi) (xi) hi \n) = j^Yl (l]|6p,5t,6k5t,^ 



k.k' 



L3 



r[(^/--Ei)ti-(p/-p.)-xi] 



and therefore: 



The third term is given by: 

T^'^HPflP'i'^t) = J '^^l J y"'^^! J dx2'W (Xi) W {X2) 

•(f}|6p^.^| {Xi)^l (Xi) (X2) ^/ (X2) 6^, |f^> . 

The matrix element becomes: 

{n\bp^4 (xi) (xi) v-l (X2) v^, (2:2) \n) = 

= ;^ ^e-*(^'=*-''-^)e^K'*^-'''"'0(f)|&p^6L^z(xi)Vj(x2)6k6Ll^) 

k,k' 

^ l^e-rJ'^^''-P'-'''^ei'^^f''-Pf'''\n\^Pi{xi)4>]{x2) \n) 

= J_ g-i(Sit2-p,-X2)gi(B/ti-p/-xi) ^g--i(£;fcti-k.xi)gi(£;fct2-k-X2)^ 



and thus we obtain: 

T^^^ (p/;p.;t) 



-7™^ dtij^'dt2j^ J d^ij dx2U.(xi)u;(x2)e-^(^'*=-P--"^)e^(^^*i-P^"^) 

, ^ g-i(Bfcti-k-xi)gi(£;fct2-k-X2)^ 



(28) 



(29) 



(30) 



(31) 



(32) 



C. Computation of the Transition Probability 

The next step is to compute the transition probability Pkj {pi]t), defined in Eq. (|16|) : now we reintroduce the 
mass-index k and j. In the interaction picture, it reads: 



Pfc,(p,;t)-^e*K' ^^'0*E[T;fe(p/;p,;t)T,,(p/;p,;t)] , 



(33) 



Pf 



where Tjj (p/;Pi;i) has been defined in Eq. ((25|) . When average over the noise, the non-zero terms are those which 
contain an even number of noises (in the graphical representation in terms of Feynman diagrams, only products of 
diagrams having an even number of dotted lines). Keeping terms only up to second order, we have: 



E[TA.(p/,p„i)r/,(p/,p,,0] = <5p,.p, +Jp„p,E[rf'(p/;p,;i)] + E [rf * (p/; p,; t) 

E rW*(p/;p,;i)T;^'^(P/;p.;t) 



^p/,Pi 



(34) 
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Therefore we get: 



E^'^'-E'^-" )t 



l + /f (p,;i) + /f *(p,;0 + 4fc' iP^■,t) 



(35) 



where we have introduced the abbreviations: 



/ " / y E 



Tr*(P/;P.;t)Tj ' (p/;p.;t) 



(36) 



p/ 



These two terms are exphcitly computed in Appendix B and C. Here we report only the final result, which is: 



2 / dsf{s){t^s) 





1 ^3/2 

{2nf r% 



(37) 



and 



if^ (pj;0 = -Imj 



dsfis) (t-s) 



1 7r3/2 



{2n) 



(38) 



where /(s) is the temporal correlation function of the noise. Inserting the expression for /■^•' {pi;t) and (jpi;t) 



written above in Eq. psp one obtains 



.(fif -_B<^"')t r _ 7 (mj - mfc) 



87r3/ 



9 S 9 



dsf (s) (t - s) 



dsf is) [t - s) 



(39) 



where in the second line we use the definition Eq. (|T3|) . Finally, the transition probability, Eq. (|16p . is computed to 
be: 



PKO^K»iP^) = \ [Pss{p^■,t)~PLs{p^■,t)~PsL{p^■:t)+PLL{p^;t)] 



1 — cos 



_ 7(mj - mfc)^ 
87r3/2r3,m§ 



ds/ (s) (t - s) 



(40) 



In analogous way we can find the probability that a K'^ remains K'^ with the time evolution: 



Pro^K" (P<) = ^ [PSS iP^■, t) + PlS {P^■, t) + PsL {p^]t) + PlL {p^]t)] 



1 + cos 



- ( e\^^ - Ef^ ) t 



_ 7(mj - rukf 



dsfis) (t-s) 



(41) 



The computations starting from the anti-kaon are completely analogous, since we neglected CV violation. Let us list 
here on two observations: 



(1) The probabilities Pk°^k° (Pi) a-nd Pxa-s^iio (Pi) sum to 1, which means that there are no particle losses. 

(2) The factor after the cosine is due to the noise present in this model, i.e. if there is no noise (7 = 0) the square 
bracket gives 1 and arrive at the standard oscillations formula^ 0, [2^ . 



In some of these works, inside the cosine, only the mass difference appears. This is the same result as ours, if one makes the approximation 
(good at the non-relativistic level) Ey' = rrijC^ + ~ rrijC^. 
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Since the effect of the noise is usually very similar to the one given by the decoherence [1, , and it is well know 
from the literature that in such cases decoherence dumps the oscillation with an exponential function, we can think 
at the term inside the square bracket as the first term of the series of an exponential function. Therefore, we can 
guess that the exact result for the probability transition from a kaon to an anti-kaon could reasonably be: 



- { eI''' - e: 



(S) 



exp ■ 



7 (mj - mfc) 
87r3/2r3m2 



dsf{s) it 



(42) 



To conclude the computation on a single particle, we can include the decay of the particle adding to the free 
Hamiltonian Hs an imaginary term: 



H. 



H. 



2^- 



(43) 



This changes the previous computation by sending E^ — > E^ — ^Tj and one can easily notice that the only change 

consists in multiplying each P^j (p,; t) with the exponential function e * in order to take the decay into account. 

Thus we obtain the final result: 



■-t , -i#t 



2 cos 



e: 



E. 



_£L±£Sy 



T-(m3-mj.) 
,fi_3/2^3 „2 



effect due to CLS model 



(44) 



where now we have assumed that the noise is white in time, i.e. /(s) = 5{s). This is the standard situation with the 
CSL model. Similar results hold for the other transition probabilities. 



IV. THE COLLAPSE MODEL FOR TWO PARTICLE STATES 

At a 0- factory, e.g. in the KLOE detector at the DA$NE collider [s^, H^-fl^, neutral kaons are copiously produced 
in an entangled antisymmetric state: 

^ \KsKl) ~ \KlKs) ^ \K^K^)-\K^k^) ^ ^^^^ 
^/2 \/2 

We want here to investigate how the mass-proportional CSL model changes the time evolution of entangled states. In 
order to obtain a more general result, we perform the computation for an arbitrary two-particle state. One particle 
evolves to the left hand side and the other kaon evolves to the right hand side, for which they need the time ti^t^, 
respectively. It is convenient to use the mass basis \KsKs) , \KsKl) , \KlKs) , \KlKl), where the state on the left 
hand side is a plane wave with momentum — and the one on the right hand side is a plane wave with momentum p^. 
In the Fock space framework, this means for example \KgK]^) = Og {—Pi) (Pi) \^)- Since we are considering only 
states with two particles with definite momentum ip;. the four states above form a complete basis. So the generic 
initial state can be decomposed as: 

^3k\K,Kk) with Yl = (46) 

j.k=S.L j.k=S.L 

Let us assume that we want to know the probability to find the left particle at time ti in the state \Fi,pi) = 
Y^m=s L^rn \KrmPi) and the right particle at time tj. in the state |i^r7Pr) = X]n=SL7n l^mPr)- Here and \Fr) 
can be mass or flavor eigenstates. In the end we have to sum over all p; and all p^ since we are interested in a result 
independent from the particular final momentum of the particles. Wc start by computing the amplitude: 

= Y {Fi,pi \Uj iti)\Kj, -p,;) (Fr, p, \Uk iU)\Kk, p,;) = 

j,k=S,L 

Y OijkP*ml*n (Km, Pi \Uj (ii)l Kj, -p,;) {K,„ Pr \Uk {tr)\ Kk, P^) = 

j,k,m.7i — S,L 

= (^jk^;i*k{Kj,Pl\U,{ti)\K„-p,) {Kk,Pr\Uk{tr)\Kk,P^) , (47) 

j,k=S,L 
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since it gives the probability of interest: 



Pi,Pr 



X E 



{Kf , Pi I Uj, [ti )\Kj,,^T^iY {Kj , Pi I Uj [ti ) I Kj , p. 

Pi 

{Kk',Pr \Uk' itr)\Kk',Pt)* {Kk,Pr Pk {tr)\Kk,Pt 



(48) 



The noise average involves terms that are contained in the left square bracket (the particle on the left hand side), 
or in the right square bracket (the particle on the right hand side) and mixed terms where the correlation is taken 
between a piece from the first square bracket and the other piece from the second one. Since we work with plane 
waves all these different terms are of importance. In the more realistic situation, however, where we assume to have 
initial confined wave packets that propagate in opposite directions, it can be show (see appendix D) that the mixed 
terms give a negligible contribution compared to the other ones. This considerably simplifies the computation to 



PiFr,Fr) = [{Kr,pi \U,, (tz)| A>,-p,)* {Kj.pi \Uj {ti)\K,,~p,)\ 

pi 

X ^ [E {Kk',Pr \Uk' {tr)\Kk',P^)* {Kk, Pr \Uk {tr)\Kk, P.)] . 



(49) 



The terms inside the square bracket are the same as those we computed in the previous sections for the single 
particle case, therefore we obtain a factorization of the probabilities, i.e. 



P {Fi;Fr) = J2 ajkP*lla*,yPy^k' Pfj {-p^; U) ■ Pk'k {P^■, U) , 
j,k,j' ,k'=S,L 



(50) 



with Pkj (pi;t) given by Eq. ((39)) . 

For the antisymmetric initial state (|45|) (asL = —ctLS ~ else ay = 0) and choosing the final states, \Fi) = j-f^") 

{Ps = Pl = and = \K^) (75 = 7l = -^), respectively, we have to compute 



1 



P{K''-K'') = -[Pssi-p^■,U)■PLL{p^■,U)+PLshP^■,U)■PsL{p^■,tr] 

+ PsL (-P,; i/) • Pls {p^■, U) + Pll {~Pi;tl) ■ Pss {Pi;tr)] 
where in the case of white noise field we have: 



and thus: 



P 



Pkj {Pi;t) = e 



(51) 



(52) 



ti-—tr _|_ g-—ti--ir*^ 



2 • cos 



(53) 



effect due to CLS model j 

In the next section, we comment on the experimental implications of the above formula. 



V. ESTIMATION OF THE EFFECT OF THE COLLAPSE ON SINGLE AND ENTANGLED KAONS & 

CONNECTIONS TO OTHER MODELS 

The effect of the collapse is the introduction of an additional damping of the interference term. This effect is also 
observed in the case one assumes that during the time evolution of the single or two-particle state the kaon or kaons 
interact with an environment in the mass basis {Ks, K^}, which was investigated by the authors of Ref. [29| . They 
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treated the neutral kaons as an open quantum system with Markovian interactions (for an overview on open quantum 
system see e.g. (soj ) and their dynamics is given by the Liouvihe-von Neumann equation with an additional Lindblad 
term: 



(54) 



which gives rise to the well known Lindbald equation [SlJ, |32 



-p = -iH p + tpH^ -V[p] 



(55) 



They further assumed that the interaction occurs in the mass basis. For the single particle case one chooses the 
Lindbald generators to be Ls^l = \/ Asingio | ^s, l ) {Ks,l \ where Agingio is the strength of the interaction. The solution 
of the components of the density matrix p = Pij{t)\Ki){Kj\ is given by: 



(56) 



thus, as we expect, only the off-diagonal terms are affected by the interaction with the environment. In particular, 
the probability to find an antikaon when a kaon state was produced at time t = Q becomes: 



e a * — 2 cos 



— {ms - mL)t 
n 



£s±£l , 



dccohcrcncc effect 



(57) 



which is of the same structure as the result, Eq. p4)) . based on the mass dependent CLS model. In the two particle 
case we have a similar behaviour if the two Lindbald generators are chosen to be Li = ^Atwo-particio|ei)(ei| and 
L2 = \J Atwo-particie 1 62) (e2 1 with |ei) ~ \K s) ® and \e2) — \Kj}i ® \Ks)- Then for the initial antisymmetric state 
p — with IV'^) — "^{Ifii) ^ |c2)} we obtain the following time dependent density matrix before one kaon 

decayed: 



P{t) 



1 „£5±£l- 

— p IK ^ 

2 ' 



ei)(ei| + |e2)(e2| 



(|ei)(e2| + |e2)(ei|) 



(58) 



Assuming that after one particle decayed no further dccoherence effect occurs, namely only the two-particle state 
interacts with the environment, then we obtain the following joint probability [29| : 



P 



2 ■ cos 



{ms - rriL) (U - U) 



(59) 



decoherenee effect 



that has the same structure as the result obtained from the CSL model, Eq. except that here the damping is not 
sensitive to the sum of the times, but on the time value corresponding to the first measured kaon. The underlying 
philosophy is the assumption that an entangled two-particle state has a single time evolution and interacts as a 
whole system with the environment. When one kaon is measured or decays there is no further interaction with the 
environment or the single particle state interacts via another interaction with the environment. Thus the longer 
an entangled pair of kaons survives, the more it interacts with the environment, namely the effect increases. So 
far, though investigations for kaons and B-mcsons (4l| have been tried, no effect has been found from the currently 
available data. 

Another option is to assume that both kaons, though being entangled, separately interact with the environment. 
This is the case for the considered CSL model, because it assumes that the contribution from the collapse is the sum 
of contributions to the Hamiltonian of the short-lived state and the long-lived state. 

Let us now compare with data of experiments and estimate the effects. There exist precision experiments for neutral 
K-mesons and neutral B-mesons in the antisymmetric entangled state \ip~) that investigate bounds on these possible 
decoherenee effects or differently stated on the stability of nonlocal correlations. They consider the joint probabilities 
and model them by one single phenomcnological parameter C multiplying the interference effect, i.e.: 



2 • cos 



■-tr 



E, 



e-^(*'+*-).(l-C) 



(60) 
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For data they get the following numbers. The authors of Ref. (SJ] analyzed the 1999-CPLEAR experiment [Sa] and 
found the following upper bounds if the decoherence effect is assumed to be in the mass basis (as discussed above) 
averaged over both measured time configurations: 



c^o.ian:^ 



(61) 



from which we estimate via 1 — C = e 



o-particic-min{t,,tr} magnitude of the averaged interaction strength to be: 



Atwo-particlo — 0.28 • 10^ S ^ 



or 



A, 



A. 



two- particle 



two-particle 



= 0.25. 



(62) 



Let us compare this result with the theoretical prediction given by the mass-proportional CSL model. Using the 
stronger value suggested by Adler for 7, namely 7 = lO^^^cm'^s"^ [23I, and recalling that rc = 10~^cm, \ms — m^l ~ 
3.5 • 10^i2MeV/c2 m, and toq - 9.4 • lO^MeV/c^ [H, we obtain: 



A 



CSL 



7 {ms - tul) 
ISTT^/^r^mg 



= 1.5 X 10 



-38^ 



or 



AcsL := ^ = 1.3-10- 



(63) 



Thus the expected effect is by many orders smaller than the sensitivity given by the CPLEAR experiment [3g . 

The KLOE collaboration |37H39| exploring the antisymmetric entangled kaon state produced at the $ resonance 
has also investigated possible decoherence effect. It is a very clean experiment to look for small effects. The above 
idea of the ^-parameter [3^ was investigated — differently from the above described experiment — for the final states 
7r+7r~ on both sides. It gives so far the highest precision due to a CV suppressing mechanism, i.e. the value is: 



= 0.003 ± O.OlSstat ± 0.006, 



'syst 



(64) 



Collapse models apply to all systems, and for small times one has: C, ~ Aqsl • t- Therefore we have to conclude again 
that the sensitivity is not high enough to show the possible effects of the CSL model. Certainly, one has in this case 
to do the calculation including CV violation, but the final conclusion cannot change significantly. 

Let us now discuss the other meson systems. For all mesons types except the K-mesons the decay widths are in 
good approximation equal (see appendix A), — — T, thus the probabilities to find on both sides the same 
fiavor eigenstates = B, D, Bs at different times is given by (compare with Eq. (|60l) ): 



P{M°,tv,M\tr) = 



e R 



1 



cos 



1 



- ( - K 



(65) 



hence the effect of the decay property is only an overall effect. Let us compute the effect predicted by the CSL model 
using the data of [s^ for these mesons: 



A 



CSL,S-mcson 



A. 



CSL-Sg-meson 



A 



CSL.D-meson -"^ 



7(TO_f/ - tul) 
167r3/2r^77iQ 

7(mg - iriLf 
167r3/2r^77iQ 

7(mg - triLf 



1.4 X 10 



1.7 X lO^^^s"^ 



3.2 X 10"^V^ 



or 



AcSL 
AcSL 



A 



CSLjS-ineson 



r 



AcSL, 



A 



A 



CSL 



CSL,D-meson 



2.1 • 10" 



-46 



= 2.6-10- 



1.3 - 10" 



-49 



(66) 



We observe that the effect of the collapse compared with the mean time is the best for the Ss-meson. For B-mesons 
estimates for (^B-mosons grists (iol - lisj , they are again far from the needed accuracy. Last but no least let us mention 
that authors discuss possible decoherence effects arising from different effects, e.g. the author of Ref. [31 discusses 
decoherence arising from the presence of dark energy in the universe and concludes that a possible decoherence effects 
due to quantum gravity strongly depends on the details of the structure of the quantum foam. 



VI. CONCLUSIONS 



Various experiments suggest that neutral mesons, particles that consist of a quark and an antiquark, are states that 
are superpositions of two different mass eigenstates. In this paper we addressed the question if models of spontaneous 
wave function collapse, which modify the standard quantum evolution by adding non-linear and stochastic terms to 
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the Schrodinger equation, are testable for single and even entangled system in High Energy Physics. Certainly, if 
spontaneous collapses occur in Nature, they would affect these systems at different energy scales, however, since they 
might be superpositions of mass eigenstates they are a unique laboratory to test collapse models. 

We have focussed our analysis on the mass- proportional CSL (Continuous Spontaneous Localization) model j2^[23j. 
In order to perform the lengthly and involved computation we had to use different methods. First of all we replaced 
the collapse noise with an imaginary noise that has been shown to produce all the experimentally testable predictions 
of the CSL model. Next, the noise has been treated as a small perturbation. Moving to the interaction picture we 
derived order by order via Dyson series the transition amplitudes. Via them we obtained the probabilities of interest. 
In the very last step we assumed that the noise is white in time. Consequently, via our computation one can in future 
also consider the effects of different noise correlation functions. 

The result of the computed probabilities within the CSL model shows that the oscillation terms are affected by 
an exponential damping. That is sensitive to the mass difference squared and the two phenomenological parameters 
of the CSL model, which, if the theory is taken seriously, acquire the status of new constants of nature. They are 
currently fixed to certain numbers [l^ to be in agreement with all current experiments. We computed also the case 
of bipartite entangled systems where we found that the effect of CSL factorizes, i.e. each meson seems to be affected 
separately and independently of the initial state they share. 

We compared the CSL model with other models [2^,[3^ that investigate possible decoherence effects or the stability 
of nonlocal correlations and that have been investigated in accelerator facilities. Our overall conclusion is that the 
possible effect on these mass-superposed states is too small to be seen in current or future accelerator facilities directly 
in the oscillations. 
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Appendix A: Kaon Phenomenology 

The phenomenology of oscillation and decay of meson-antimeson systems can be described by nonrelativistic quan- 
tum mechanics effectively, because the dynamics is depending on the observable hadrons rather than on the more 
fundamental quarks. A quantum field theoretical calculation showing negligible corrections can e.g. be found in 
Refs. 0,111. 

A neutral meson Mq is a bound state of quark and antiquark. As numerous experiments have revealed the particle 
state Mq and the antiparticlc state Mq can decay into the same final states, thus the system has to be handled as a 
two state system similar to spin i systems. In addition to being a decaying system these massive particles show the 
phenomenon of flavor oscillation, i.e. an oscillation between matter and antimatter occurs. If e.g. a neutral meson is 
produced at time t = the probability to flnd an antimeson at a later time is nonzero. 

The most general time evolution for the two state system M° — M° including all its decays is given by an inflnite- 
dimensional vector in Hilbert space: 

\m) = + m\M') + c{t)\h) + d{t)\h) + ... (67) 

where /, denote all decay products and the state \ip{t)) is a solution of the Schrodinger equation {h = 1): 

^IV-W) = -^H\m) (68) 

where H is an infinite-dimensional Hamiltonian operator. There is no method known how to solve this infinite set 
of coupled differential equations affected by strong dynamics. The usual procedure is based on restricting to the 
time evolution of the components of the flavour eigenstates, a{t) and b{t). Then one uses the Wigner-Weisskopf 
approximation and can write down an effective Schrodinger equation: 

^Mt)) = -iHm)) (69) 
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where is a two dimensional state vector and H is a non-hermitian Hamiltonian. Any non-hermitian Hamiltonian 
can be written as a sum of two hermitian operators M,T, i.e. H = M + ^T, where M is the mass-operator, covering 
the unitary part of the evolution and the operator F describes the decay property. The eigenvectors and eigenvalues 
of the effective Schrodinger equation, we denote by: 

H \M,) = A, \M,) (70) 

with Xi = rrii + ^Ti (c = 1). For neutral kaons the first solution (with the lower mass) is denoted by Ks, the short 
lived state, and the second eigenvector by Kj^, the long lived state, as there is a huge difference between the two 
decay widths Ts — 600Fi. For B-mesons the lower mass solution is denoted by with L for light, and the second 
solution by Bh with H for heavy. For this meson type (as for all the other ones except K-mesons) the decay widths 
are in good approximation equal, i.e. F^ ~ Th- Thus, the huge difference in two decay widths is special to K- mesons 
and this is one reason that they are attractive to various foundational tests, such as e.g. tests for nonlocality [46l - l5ll | 
or the very working of a quantum eraser [s^, Is^ - IH^ ] . 

Certainly, the state vector is not normalized for times t > due to the non- hermitian part of the dynamics. 
Different strategies have been developed to cope with that. In Ref. [s^ the authors developed based on the open 
quantum formalism a framework that show that the effect of decay is a kind of decoherence. In this paper we handle 
only the surviving part of the evolution, but given the framework developed in Ref. [ssj (or similar in Ref. [sij ) one 
can straightforwardly obtain the full quantum information content of the meson systems, e.g. to study Heisenberg's 
uncertainty in relation with CV violation in the time evolution [stI . [ssj . 



APPENDIX B: Computation of ij}^ {pi;t) 



Let us focus our attention on the stochastic average, i.e. 



ri^)*(p/;p.;OrW (p/;p.;i) 



E 



-_Lyr',,^|,,^,„,,^„-*l(.S"-.,")..-.P,-......l) 



y^lm,lm.kj^ / 1^*1 / j dxi j dXjE [w (xi 



(71) 



The average over the noise is: 



E[W {XI)W{X2)] =f{tl-t2) 



-(xi-X2)V4rg 



3 ' 



(72) 



where / (ti — ^2) is a generic correlation function characterizing the noise. Therefore we have: 



E 



T«*(p/;p,;i)T|'^ (p/;p,;t) 



(1) 



Jo 



dt2e " 



f{tl-t2) 



=C{t.pf,pO 
If f g-(xi-X2)V4'-c , 

— / dxi / dx2 ^ ^ ^3 e-«(PJ-P')-(^^-^^) 



47rrc) 



(73) 



= S(p/,pi) 



We first compute the function S: 



1 /■+- /•+- e 
S'(p/,p,) = -^ / d-y-i I dy.2 



-(xi-X2)V4r2 



■i(P/-Pi)-(xi-X2) 



47rrc) " 



(74) 
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Here we wrote explicitly the integration extremes (given by the box of side L). We make the following change of 
variables: 



y = (xi + X2) and x = (xi — X2) . 
The Jacobian of this transformation is 1/2^ and the integral (in the one dimensional case) changes as follows: 



dxi / dx2f{xi,X2) 



dx 



-(L-x) 



dy [/ {x,y) + f {~x,y)]. 



-{L-x) 



In our case we have: 

/(x,y) = -^- 

(V47rrcj 

Accordingly, the function S becomes: 



^(P^-P-)-"^/(x)+/(-x) = 



47rrc) 



r2 cos 



- (p/ - p.) • X 



(75) 



(76) 



(77) 



S (P/,P») 



1 1 



dx 



dy- 



r2 COS 



(P/ - pO • X 



^ / dx- 



t2 cos 



^ (P/ -Pj) -x 



(L - x^) 



23-11- i 

4=1 



(78) 



In the limit L ^ go (which we will do later) the term ^ Ili^i 2 



1 rr3 o [L-xj) 



S (p/,Pi) 



r2 COS 



(P/ - P») • X 



1 

1? 



1, so we can write: 



-L -xV4r^ 

dx-^ ^P,7r(p/-P.)->'_ 



-i (^/4^rc)' 



(79) 



We now compute the function C(i,p/,Pi). Let us introduce the two variables a = yE^p — e\''^ ] e b 
— i — e['^^\ and so we can write 



C(t,P/,P,) - f dh f dhe^'^'^e^'-fih-h) 

Jo JQ 

' dh f dt2e^[('^"^)(*i-*^)+(''+'''(*i+*=)l/ (ti - t2) 



(80) 



We now make the change of variables introduced when computing 5; we define m = ii + ^2 and s = ti — t2 and focus 
our attention only on the case b = —a: 



C(t,p/,p,) = f dh f dt2e'»(*i-*=)/(<i-<2) 

jQ JO 



2t+s 



dse^^'fis) / du+ dse'^'fis) 



2t-s 



du 



dse""-" f (s) (t + s) + dse""' f {s) {t ^ s) ^ 2 I ds cos (as) f (s) (t - s) 
Jo Jo 

In the white noise limit / (s) = S (s) and the function C takes the very simple expression: 

Cn, (i,P/,Pt) = t . 

Collecting all pieces together, we have: 



E 



T^i^'* (P/iPiii)?")'^ (P/;Pj;i) = V7m,7mfcC'(i,p/,Pj) / dx 



.(1) 



1 



L3 



>^(P/~Pi)-x 



(81) 



(82) 



(83) 
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and therefore 1^^} {pi',t) becomes: 



Pf 



dx- 



e -xV44 



47rrc)" 



ei(p/-Pi)-x. 



(84) 



We now arc in the position to take the hmit L — > oo. so that the above equation becomes: 



1 



{2Trh:f 

The last integral can be computed exactly: 



4-00 



dx- 



g^(P/-Pi)-x 



C(i,P/,Pi) 



(85) 



+ 00 



-p h 



(P/-Pi)-x ^ g- 



therefore we have: 



(27r;i)' 



(86) 



(87) 



Now, the Gaussian function has a spread <t ^ h/rc — 12 eV/c, which is very small compared to the typical 
wavelengths entering the oscillatory terms. Indeed, the module of the exponent of that term is: 



P? 



2toi. 2m 



4^^^^^P^-2,7xlO--(eV/c)-. 
2n mjurik 



(88) 



In the case j = k this is simply zero while for j ^ k we substitute \ms — — 3.5 • 10~"'^^MeV/c^, ~ ms — 
498 MeV/c^ [IHl and, for the kaons at DA$NE we have computed t using the maximum possible length L = 10 m 
and taking w = 0, 2c so that t = L/w = 1.6 • lO^^s. Regarding the term C(f,p/,Pi) it can be shown, starting by 
its definition, that it is composed of terms that involve exponential oscillating with a phase of the same order of 
the one above multiply with the Fourier transform of the noise correlation function / (t). This means that, like the 
exponential case we studied just above, even this term changes little the Gaussian width and so we can bring all these 
terms outside the integral and set p/ = Pi. In this way, we need only computing C (t, p;, pi), which is the expression 
as given in Eq. (|81|) with a = —b = 0. So it remains only a Gaussian integral that can be compute explicitly and we 
get the final result: 



dsf (s) {t - s) 



1 TT 



3/2 



J {2n) 



In the white noise case f (s) — 6 (s), which we are mainly interested in, the above formula reduces to: 



3 ^3 



(27r) 

We can notice that this formula is consistent with the relativistic expression first derived in Ref. [5£ 



.nijinkC 



■4 t 7r3/2 



"^iJ^ijf) (2.) 



(89) 



(90) 



(91) 



since, in the non relativistic limit: £'|-'-' = ^/p^c^ + mjc' 
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APPENDIX C: Computation of /f {pi;t) 



We now compute /j^'' (p^; t) = E T^^"^' (p^; p^; t) , with Tj''' (p^; Pi; t) as given by Eq. ([32]). Using once more relation 
Eq. dill), we obtain: 



.(2) 



-Ira, J2 f dti f 
1^ Jo JO 



Jo 



C"(t,Pi,k) 

g-i:(pi-k)-(xi-x2 



(92) 



=S(p.,k) 



The function 5 is the same as the one defined and computed in the previous section. The function C" instead is 
sUghtly different from the function C previously introduced, the main difference being that in C the second integral 
goes from to ti, while for C it goes from to t. As before it will be interested just in the case pi = k and, 
performing similar computation, one finds out that: 



C'it,p„p,)= f dh f dt20 {h ~ t2) f ih - h) = f ds/(s)(t-s) 
Jo Jo Jo 



Going back to the definition of P- ' (pi;t), by taking the limit L — > oo, one obtains: 



/f^(p,;0 =-7m, / dkC'(t,p,;,k) 



(93) 



(94) 



where the Gaussian integral has already been computed in the previous Appendix B. At this point, as we did in the 
previous section, we can make the approximation, inside the integral, C {t, Pi, k) ~ C" (t, p;, pi) and so we get: 



if' (Pj,Si;i) = -7m,- 



dsf{s) {t~s) 



1 7r3/2 
(27r)^ rl 



whit 



c noise case 



7m, t 7r3/2 
' 2 {2i,f rl ■' 



that is still consistent with the relativistic result: 



' ^P-^-*^- 2^'"^- ^20-) (2^)3 ,3 



(95) 



(96) 



APPENDIX D: A computation with wave packets 

Let us start from the formula Ec^. (|^5|) : 

P(,Pr j,kj' ,k'=S,L 



X E 



E P' I^J' (^01 -P»>* {Kj,Pi \Uj iti)\ Kj,-v, 
Pi 

^{Kk',Tpr\Uk' {tr)\Kk',V,)* {Kk,Vr\Uk{tr)\Kk,Vt 



(97) 



Here we want to make clear why we can neglect all the terms involving the correlation between terms in the first 
square bracket with terms in the second square bracket. We recall that: 



(i^„p/|[/,(t)|i^„p,) = 



= e " / 



'5p/,P. +7;''' (p/;p.;t) + Tf' (p/;p.;t) 



(98) 
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If we substitute this in the second and third hnes and we take only the contribution up to the second order we get: 



E. 



^ (A'fc',Pr \Uk' {tr)\ Kk',Pi)* {Kk,Pr Pk (tr)| Kk,Pi) 



L Pr 



• ^ " {l + E [rj,'^* (-P,; -p,; to] + E [rf ^ (-p,; -p,; i/) 



E T*r' (p.;p.;t.) +E T^'^ (p,; p,; i. 



-E 



■E 



(1) 



+E 



tI'^ i~P^;^P^■,tl)T|^}^*iPf,p,■U) 



pi 



tIP* {pl;^p^;tl)Tl'^ {pl;-p^■,tl) 



(99) 



We can see that the first and the last two lines contain pieces that involve only terms referring to the same particle (left 
or right). But there are also pieces containing the product of terms which refer to one particle and terms which refer to 
the other particle. This pieces seems unphysical, especially because they are present even if we take a separable state 
as the initial state. And the fact that the correlation between left terms and right terms could play an important role 
even for separable states sounds suspicious. Indeed this problem arises because we had oversimplified our analysis: 
by using plane waves as states of the particles, we are working with totally delocalized states, which is the origin of 
the the correlation between terms of the left and the right particle. Anyway this is not what actually happens in 
laboratories: in such cases, the state of the system is described by wave packets well localized in space and we will 
show now that, with such an assumption, the suspicious pieces cancel. To do that, we have to replace this plane waves 
with wave packets. Therefore, in place of the two initial states \pi) and |— P;), we now consider: 



\f) =^f iP^)\P^) and Iff) = ^g(Pz) |p») , 



(100) 



where |/) {\g)) it is a wave packet propagating in the left (right) direction. Regarding the final states, we will continue 
to take them as plane waves, because, as before, we will sum over them. Taken different initial states implies that the 
matrix element {Kj,pf \Uj {t) \ Kj,pi) become: 



{K,,Pf \U, {t)\ K,J) = ^ / (p.) 



"Pf,Pi 



-T '> (p;;p,;t) + r ^^ (p/;p.;t) 



1 noise 



2 noise 



(101) 



and the same happen with the matrix element containing \g). Let's see how change the bad terms, for example the 
one that before was: 



T^'H-P^■,-P^■,tl)Tt'>* {p^■,P^■,tr) 



(1)* 



(102) 



18 



Looking at the formula for the correlation Eq. ((99)) . we know that this term came out taking the delta (i.e. the 
zero order contribution) for the bracket labeled by f and k and taking the first order term T^^") for the ones labeled 
by the index j and k' . If we take now the corresponding terms we get /* (p;), g (Pr) instead of (5p,.pi, i^p^.p; and 
Epi /(pO^i^^ (Pi;Pi;^i). Ep2.9* {P2)tI:}^* (Pr;P2;ir) instead of Tj^' (p,;pi;t;), Tj^}^* {p,;p2]tr) and so it becomes: 



E}- '-E, 



".r(p/)5(Pr) E /(Pi).9*(P2)IE[rf^(pz;pi;tz)4^'*(p,;p2;i.)J ■ (103) 

Pi.Pr Pl,P2 

Regarding this expression the interesting piece is: 



z=Y, ./(Pi).9*(P2)iE[TW(Pi;Pi;i04'^*(Pr;P2;ir] 



(104) 



Pl>P2 



dti / dxiw (xi) {Q\ '(/'!/ i^i) i^ji i^i) I/) X 



Using Eq. 

X y dt2 j dyi2w[x2){Vl\bk:p^ii\,j[x2)iik'i[x2)\g) 
We can focus on compute the matrix elements. For example the second one is: 

{^\ bk'p^^lj {X2) ^k'l {X2) \g) - (^^1 bk'p,.4'i e^^«*^ Vfe' (X2) e-^"^'^ \g) 



(105) 



(106) 



Now we can introduce \g,t2) = e"'^-^^*^ ^jjg^^ jg j^g^ wave packet evaluated up to the time ^2 with free 
evolution. From the experiment we know that the kaons arc always found one at left one at right. This means that 
the wavcpackets associated to these kaons docs not spread too much compare to their distance from the source, i.e. 
gix,t2) = (x 15,^2)-'' 

Now we can use the key relation: 



\g,t2) = / dxg{x,t2)i'l,{x)\n) 



and so: 



(107) 



{n\ bk'p^i^lj i^k'i {X2) \g) = / rfxg (x, t2) (0| bk'p^i^lj e*^^*^ Vfc' (x2) 4' (x) m = 



g{x2,t2){n\bk'p^iij{x2)ei"''^ = .g(x2,t2) 



j-(Ei'''\2-Pr-^2 



Doing the same for the other matrix element and inserting everything in Z wc get: 



= ^/Imjlm^, J dti J dt2 J dxi J dx2 

X / (xi,ii)g (x2,t2)E [w (.Ti) w (^2)] 



^^(B,"'t,-p,.x,)^i(B(''')t.-p..x: 

1? 



and using the correlation (in the white noise case) E [w (xi) w (2^2)] — <^ (^i ^ ^2) 



-(xi-X2)^/4t- 



(108) 



(109) 



One can claim that the wavepacket observed in the laboratory is the one given by the complete evolution, taken into account also the 
noise, and so it is different from g (x, t) that is the wavepacket evolved with the free Hamiltonian. Anyway it is well know from the 
literature that the effect of the noise on the spread of the wavefunction in typical times of this experiment is negligible, so we can be 
sure that if the real wavepacket remains confined, the same happens also for g(x,t). 



19 



mini 



X 



dxi / dx2e-^P' ''ie-^P'-''V(xi,ii)5(x2,ii) — —■ (HO) 

We recall that rc = 10~^m so, a part for the starting instants, the two wave packet are always far away and the 
integral is small. This doesn't happen for the terms where the correlation is taken between noise related to the same 
particle, because in such case we end up with terms containing pieces hke / (xi, fi) / (x2, ti) or g {xi,ti) g {x2,ti) 
inside the integrals, that gives an important contribution even for large time. 
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